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Abstract 


The spectroscopic constants and vibrational frequencies for the 1 3 A states of HNC, 
DNC, HOC + , and DOC + are computed and discussed in this work. The reliable CcCR 
quartic force held based on high-level coupled cluster ab initio quantum chemical com- 
putations is exclusively utilized to provide the anharmonic potential. Then, second- 
order vibrational perturbation theory and vibrational configuration interaction meth- 
ods are employed to treat the nuclear Schrodinger equation. Second-order perturbation 
theory is also employed to provide spectroscopic data for all molecules examined. The 
relationship between these molecules and the corresponding 1 3 A HCN and HCO + 
isomers is further developed here. These data are applicable to laboratory studies in- 
volving formation of HNC and HOC + as well as astronomical observations of chemically 
active astrophysical environments. 

Keywords: Astrochemistry; Quartic Force Fields; Coupled Cluster Theory; Vibra- 
tional Perturbation Theory; Vibrational Configuration Interaction Theory; Vibrational 
Frequencies; Triplet States 


Introduction 

The presence of HNC and the isoelectronic HOC + has been known in the interstellar medium 
(ISM) for some time. 1-5 HNC has been hypothesized as a product in the dissociative recombi- 
nation of interstellar 6 HCNH + in molecular clouds 7,8 as one possible astronomical synthetic 
pathway. High-resolution experimental studies have recently examined the vibrational fre- 
quencies of HNC and the deuterated isotopologue, DNC. 9-11 Additionally, highly-excited 
vibrational lines of HNC were observed toward the carbon-rich star IRC +10 216 with the 
Herschel Space Observatory 12 before the space-based telescope closed its eyes in the spring 
of 2013. 

There has been speculation that the lowest energy triplet states of HCN, HNC, HCO + , 
and HOC + contribute to rotational signals observed in the THz spectral range of labora- 
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tory studies designed to simulate various astronomical environments. 13-16 Previous CASPT2 
computations report a 4.47 eV 1 3 A' <r- X 3 S + excitation energy for HNC (Ref. 17), which is 
similar in energy for the same excitation in HCN at 4.69 eV (CASPT2/ANO-L from Ref. 18). 
Additionally, the energy separation between 1 3 A' HCN and 1 3 A' HNC has been computed 
to be 10.9 kcal/mol [CCSD(T)/aug-cc-pVTZ] and even less between 1 3 A' HCO + and 1 3 A' 
HOC + at 2.5 kcal/mol. 19 As a result, overlap in the rotational signals for the two isomers of 
the respective nitrogen- or oxygen-containing species generated is likely. Resolution of their 
individual spectra would benefit from reference data in order to distinguish the spectral 
features. 

Quartic force field (QFF) computations of vibrational frequencies and spectroscopic con- 
stants for the X 1 E + states of both HNC and HOC + have long served as examples of the 
accuracy that quantum chemical calculations can provide. 20,21 Building from that framework 
and the two decades of development since then, 22-26 the largely unexamined vibrational fre- 
quencies and spectroscopic constants for the 1 3 A' states of HNC and HOC + are investigated 
here, whereas earlier rovibrational studies have produced reference data for the 1 3 A' states 
of HCN and HCO + . 19 Additionally, with continuing growth in infrared (IR) and far 1R 
telescopic capability, reference data for these molecules could assist in the chemical charac- 
terization of astronomical spectra taken with instruments on plane- or space-based telescopes 
such as the Stratospheric Observatory For Infrared Astronomy (SOFIA) and the upcoming 
James Webb Space Telescope (JWST), as well as the ground-based Atacama Large Millime- 
ter Array (ALMA). 


Computational Details 

The procedure for the computation of the vibrational frequencies and spectroscopic constants 
of 1 3 A' HNC and HOC" 1 " closely mirrors that used for the study of the 1 3 A' states of HCN 
and HCO + (Ref. 19) and used in the study of other systems. 22,23,27-29 Due to the open- 
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are 


shell nature of these triplet systems, restricted open-shell reference wavefunctions 30-32 
combined with the coupled cluster singles, doubles, and perturbative triples [CCSD(T)] 
level of theory. 33 The geometries are optimized with Dunning’s correlation consistent aug- 
cc-pV5Z basis set 34-36 refined by corrections for core-correlation from the Martin- Taylor 
(MT) basis sets designed to properly treat inclusion of the core orbitals in the CCSD(T) 
computations. 37 From these reference geometries, central difference displacements of 0.005 
A for the bond lengths and 0.005 radians for the bond angles are generated from the simple- 
internal coordinates leading to a grid composed of 129 points. Since this is a fourth-order 
Taylor series approximation to the anharmonic potential function, the displacements are 
treated up to fourth-order to give the QFF. At each point on the QFF surface, a three- 
point complete basis set (CBS) limit extrapolation 38 energy is computed from utilization of 
Dunning’s augmented triple, quadruple, and quintuple zeta basis sets. To this CBS energy, 
corrections for core-correlation and scalar relativity 39 with the cc-pVTZ-DK basis set are 
included to form the CcCR (CBS ■$ core-Correlation + Relativistic) QFF. All electronic 
structure computations employ the MOLPRO 2010.1 quantum chemistry package. 40 

The equilibrium geometry and the force constants are produced from a least-squares fit- 
ting of the points in the QFF where the sum of squared residuals is on the order of 10 -1 ' 
a.u. 2 . Cartesian derivatives are subsequently determined using the INTDER program, 41 
and fed into the SPECTRO program 42 where vibrational second-order perturbation the- 
ory 43-45 returns the spectroscopic constants and vibrational frequencies. Different from the 
second-order vibrational perturbation theory (VPT2) results, the vibrational configuration 
interaction (VCI) frequencies are produced from the MULTIMODE program 46,47 after the 
QFFs are transformed into a Morse-cosine coordinate system. 25,48 
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Results and Discussion 


Geometrical Data and Spectroscopic Constants 

The difference in CCSD(T)/aug-cc-pV5Z energy for the optimized geometries of the X 1 £ + 
state and the 1 3 A' state of HNC gives an adiabatic excitation energy of 4.64 eV, similar 
to the 4.47 eV transition energy computed previously with CASPT2. 1 ' This is a n — > n* 
excitation out of the highest energy n bonding orbital in linear HNC into the lowest energy 
7 r antibonding orbital. As such, a (core)3a /2 4a /2 5a /2 la //2 6a /1 7a /1 configuration is produced in 
the bent triplet from the (core)3o‘ 2 4(T 2 5o' 2 l7r 4 27r 0 orbital configuration present in the linear 
singlet. The same process occurs in the HCN isomer, as well, with a similar CCSD(T)/aug- 
cc-pV5Z adiabatic excitation energy of 4.84 eV. 19 Previous adiabatic computations have 
placed the energy difference between the 1 3 A' states of HCN and HNC between 10.4 and 
10.9 kcal/mol. 19 Here, we refine this non-zero-point value to 10.3 kcal/mol from the more 
complete CcCR energies of both triplet isomers. 

The 1 3 A' HNC equilibrium geometry fit from the CcCR QFF given in Table 1 is in line 
with that produced from previous CASPT2 computations 17 but differs enough to show the 
refinement brought about by the higher-order computations employed here. The dipole mo- 
ments have a larger variance between the two methods, 1.06 D from CCSD(T)/aug-cc-pV5Z 
computations based on the CcCR geometry and 1.2963 D from the CASPT2 computation 
and geometry, but they differ by as much as one would expect between second-order pertur- 
bation theory and coupled cluster quantum chemical methods. For comparison, a previous 
CCSD(T)/aug-cc-pVTZ optimized geometry dipole moment of 1 3 A' HNC is 1.05 D. 19 When 
analyzed with respect to the CcCR 1 3 A' HCN structure, HNC has an understandably dif- 
ferent geometry, especially for the bond angle which decreases from 120.587° for ZH— C— N 
to 112.151° for ZH— N— C. Zero-point (R a ) vibrational averaging for the HNC geometry 
increases each of the bond lengths as well as the bond angle as compared to the HNC equi- 
librium geometry. The B- and C-type rotational constants are subsequently decreased upon 
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vibrational averaging whereas the H-type constant increases from 624.286 4 GHz to 626.373 
7 GHz when the vibrational effects are included. 

Since HOC" 1 " is isoelectronic to HNC, the 1 3 H / state results from the same excitation 
scheme as that given above. The CCSD(T)/aug-cc-pV5Z adiabatic 1 3 ZL' •(— X 1 X + exci- 
tation energy for HOC + is 4.21 eV. This is substantially less than the 5.86 eV adiabatic 
excitation energy for HCO + computed in the same manner. 19 The CcCR energy difference 
between the 1 3 ZL' states of HOC" 1 " and HCO + is a mere 1.4 kcal/mol, 0.2 kcal/mol less than 
CCSD(T)/aug-cc-pV5Z computations at the optimized CCSD(T)/aug-cc-pVTZ geometries 
previously indicated. 19 However, the barriers to isomerization computed previously 19 indi- 
cate that rearrangement of the isomers has a significant energy cost at more than 40 kcal/mol 
for both HCO+/HOC+ and HCN/HNC. 

The 1.80 D dipole moment for 1 3 H / HOC" 1 " (also given in Table 1) is smaller than the 
corresponding HCO + isomer at 2.51 D but nearly coincident with the CCSD(T)/aug-cc- 
pVTZ 1.81 D dipole computed as part of our earlier study. 19 The B- and C-type rotational 
constants (whether from equilibrium or vibrationally averaged) for 1 3 ZL' HNC and HOC" 1 " 
are equivalent between the molecules to within 0.5 GHz. However, the ZL-type constants 
should allow for them to be spectrally distinguished. Deuteration of either molecule studied 
here noticeably changes the rotational constants, as one would expect. 

The symmetry- internal quadratic, cubic, and quartic force constants are given in Table 2. 
The gradients are zero by construction. The force constants given in Table 2 are numbered 
sequentially from the geometrical parameters given in Table 1 where coordinate 1 corresponds 
to the X— C bond, coordinate 2 corresponds to the X— H bond, and coordinate 3 corresponds 
to ZH— X— C with X = N or O. Table 3 contains the other spectroscopic constants including 
the vibration-rotation interaction constants, the quartic and sextic centrifugal distortion 
constants, and the Watson ^-reduced Hamiltonian terms for both 1 3 ZL' HNC and HOC" 1 " as 
well as their respective deuterated isotopologues. 
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Vibrational Frequencies 

The fundamental harmonic and anharmonic vibrational frequencies as well as the first over- 
tones and combination bands of 1 3 A' HNC, DNC, HOC + , and DOC + are listed in Table 4. 
The VPT2 1 3 /V HNC and HOC + computations require input of a 2z/ 2 = u \ Fermi resonance. 
VPT2 for both DNC and DOC + requires input of a rq + z/ 2 = zq Fermi resonance. No Cori- 
olis or Darling-Dennison resonances are present in these systems. The MULTIMODE VCI 
computations are complete and exact at the three mode representation (3MR) level since all 
fundamental vibrational modes for these non-linear triatomic molecules can be coupled at 
this level. Basis set convergence is defined here as a change of less than 1.0 cm -1 with an in- 
crease in more than 100 basis functions. Convergence is achieved for each of the molecules in 
this study upon the inclusion of 420 basis functions in the computation. The bond stretching 
bases are each composed of 36 Gaussian integration points, 31 primitive harmonic oscilla- 
tor basis functions, and 20 Hermite-Gauss (HEG) quadrature points. The bending bases 
contain 26 Gaussian integration points, 23 primitive harmonic oscillator bases, and 14 HEG 
points. These basis setups are of the same size or larger as those employed to study larger 
tetraatomics where the hydrogen stretching frequencies, for instance, have been computed 
to be within 5 cm -1 or even better than 1 cm -1 as compared to gas phase experiment. 23,49,50 
The CcCR QFF is the only QFF used in this study since previous work has shown this QFF 
to be the most cost-effective means of computing highly-accurate rovibrational reference 
data. 19,22,23,27,28 

The CcCR 1 3 H' HNC harmonic frequencies from Table 4 highlight the difference in 
geometry brought about from isomerization originating with 1 3 H' HCN. N— H stretches 
occur at higher frequencies relative to C— H stretches, and such is highlighted here. Even 
the C— N bond stretch varies between the isomers decreasing from 1589.5 cm -1 in 1 3 H' 
HCN to 1568.9 cm -1 in 1 3 H' HNC. Previous CASPT2 computations 17 for 1 3 H' HNC report 
a lower harmonic frequency N— H stretch at 3103 cm -1 relative to the CcCR 3189.3 cm -1 
harmonic frequency. The z / 2 and za> modes have higher harmonic frequencies in the previous 
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work 17 compared to the CcCR QFF. 

The VPT2 and VCI anharmonic frequencies should be more accurate than either of the 
sets of harmonic frequencies reported. The zq fundamental and the combination bands and 
overtones containing it are the most anharmonic of the frequencies by magnitude. VPT2 
and VCI treat the nuclear Schrodinger equation through different means, but they arrive at 
similar if not identical frequencies for 1 3 A ' HNC when using the CcCR QFF. The N— C zq 
stretch is computed to be 1534.0 cm -1 for both methods. Further consistency is present in 
the hydride stretch where VPT2 reports this mode to be 2939.9 cm -1 while VCI is exactly 

1.0 cm -1 less. The bending mode has a larger difference between the two methods, but it 
is still only 2.2 cm -1 . Larger variances between the methods begin to arise in the overtones 
and combination bands. This is epitomized in the zq + z/ 2 combination band. VPT2 reports 
a 4456.2 cm -1 frequency while VCI reports this frequency to be 4447.1 cm -1 . However, this 

9.1 cm -1 difference is fairly small since combination bands are more difficult to model with 
QFFs than the fundamentals. Regardless, the consistency in the methods indicates a marked 
measure of reliability for these frequencies. 

1 3 A' DNC continues the consistency reported for the standard isotopologue. Within the 
Born-Oppenheimer approximation used in the formulation of the QFF, the force constants 
are the same between 1 3 A' HNC and DNC. VPT2 and VCI consider the different effects of 
the heavier atom in solving the nuclear problem. As such, each of the fundamental vibrational 
frequencies for DNC are less than their corresponding HNC frequencies. The z/ 2 stretching 
frequency is identical between VPT2 and VCI even though it is 8.5 cm -1 lower in frequency 
at 1523.5 cm -1 after deuteration. The zq and z/ 3 modes are actually more consistent between 
VPT2 and VCI for DNC than they are for HNC. Additionally, VPT2 and VCI are in much 
better agreement for the overtones and combination bands of DNC, as well. The largest 
difference between the methods is 5.2 cm -1 , again for the zq + z/ 2 combination band, but 2zq, 
2 z/ 2 , and z/ 2 + zq have VPT2-VCI differences of 1.5 cm -1 or less. 

The zq O— H stretch of 1 3 /V HOC + is the highest frequency mode of all the fundamentals 



reported in this set of triplet triatomics. It is 3146.6 cm -1 from CcCR VPT2 and 3146.0 
cm -1 from CcCR VCI, a further indication that the frequencies are well described with either 
method. The u 2 O— H stretch differs by 0.1 cm -1 here, and the z/3 bending frequency differs 
by 1.6 cm -1 with VCI producing a 1166.8 cm -1 frequency for this mode. Again, the CcCR 
QFF is providing consistent results with either vibrational calculation method. The 2zq and 
2z/ 2 overtones are quite consistent between the methods as they were for DNC. The largest 
VPT2-VCI difference for HOC + for any of the fundamentals, overtones, or combination 
bands is the 2rq bending overtone. VPT2 gives 2306.1 cm -1 while VCI is 2299.5 cm -1 , a 
difference of 6.6 cm -1 . 

Deuteration of 1 3 A' HOC + to yield 1 3 A' DOC + also reduces the vibrational frequencies 
as compared to the standard isotopologue. All of the VPT2 fundamental frequencies, how- 
ever, are within 1.0 cm -1 of their VCI counterparts. The zq + z/ 3 combination band represents 
the largest VPT2-VCI discrepancy with the frequencies differing between the two methods 
by 6.3 cm -1 . However, the v 2 fundamental and 2z/ 2 overtone are computed to possess nearly 
the same frequencies wether VPT2 or VCI is employed. 

Conclusions 

Since the X 1 £ + forms of HNC and HOC + are known to exist in the ISM, and it is rea- 
sonable to believe that the 1 3 A' forms of HNC and HOC + can arise in various interstellar 
environments, an understanding of their vibrational and rotational properties would serve to 
better inform astronomical observations and laboratory studies where they may be found. 
The corresponding isomers of each (i.e. , HCN and HCO + ) have been examined previously, 19 
and this study represents a completion of the rovibrational analysis for minima on each of 
the respective potential energy surfaces. The 1 3 A' HCN and HCO + isomers are shown here 
with highly-accurate methods to be only 10.3 kcal/mol and 1.4 kcal/mol, respectively, lower 
in energy than the corresponding HNC and HOC + isomers, though it was shown previously 
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that there is a large barrier (> 40 kcal/mol) separating the respective minima. Thus, if 
triplet HNC or HOC + are formed, they will not easily isomerize to the lower energy HCN 
or HCO + isomer. Indeed, there has already been speculation that the lowest triplet states 
of HCN, HNC, HCO + , and HOC + might be present in laboratory experiments designed to 
simulate astrophysical conditions. The 1 3 A' X 1 £ + HNC excitation is 4.64 eV and lower 
at 4.21 eV for HOC + . However, if 1 3 A' HNC and HOC + (as well as HCN and HCO + ) 
result from a formation mechanism favoring triplet states in astronomical regions where the 
molecular density and collisional probability are low, then they could exist for a long time 
lacking a radiative mechanism to the ground electronic singlet state. 

The fundamental vibrational frequencies for both 1 3 A' HNC and HOC + , as well as 
the deuterated isotopologues are reported in this study. The set of frequencies for each 
molecule are in very good agreement between the VPT2 and VCI methods with both making 
use of the CcCR QFF. Additionally, the overtones and combination bands are also fairly 
consistent giving a strong indication that the computed results presented here are reliable. 
The spectroscopic constants provided from second-order perturbation theory also give the 
potential for further insights into the behavior of these molecules. All could be of benefit for 
analysis of spectra whether taken in the laboratory or from an astronomical observation. 
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Table 1: The CcCR QFF 1 3 ZL' HNC and 1 3 ZL' HOC + Minimum Energy Structures 
(in A and degrees) and Rotational Constants (in GHz). 


Molecule 

Equilibrium 
This work Previous" 

Zero-Point 

Standard Deuterated 

1 3 ZL' HNC 

R(N-C) 

1.261 56 

1.2709 

1.265 80 

1.266 15 


R(N-H) 

1.032 80 

1.0417 

1.049 00 

1.045 17 


ZH-N-C 

112.151 

111.63 

112.325 

112.235 


Zf e /0 

624.286 4 


626.373 7 

356.142 0 


B e / o 

44.809 22 


44.696 35 

40.840 85 


Ce/O 

41.805 09 


41.504 93 

36.419 59 


h 

1.06 

1.2963 



1 3 ZL' HOC+ 

R(O-C) 

1.224 90 


1.229 27 

1.229 50 


R(O-H) 

0.993 92 


1.008 66 

1.005 15 


ZH-O-C 

115.780 


116.119 

115.959 


Zf e /0 

708.961 6 


717.653 0 

404.429 5 


B e /o 

44.945 07 


44.782 93 

41.085 23 


C e / o 

42.262 07 


41.943 11 

37.077 49 


/' 

1.80 





“10/(8, 2) CASPT2 computations from Rajendra and Chandra (Ref. 17). 
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Table 2: The CcCR QFF Quadratic, Cubic, and Quartic Force Constants (in 
mdyn / A n rad m ) for the simple-internal coordinates of 1 3 A ' HNC and HOC + . 


HNC 

HOC + 

F u 

9.641 120 

Fun 

369.16 

7*n 

11.387 673 

7*im 

431.05 

7*21 

0.383 265 

F2111 

22.85 

7*21 

0.331 626 

7*2111 

19.34 

7*22 

5.656 329 

F2211 

- 16.99 

7*22 

6.360 166 

7*2211 

- 10.40 

F:n 

0.638 345 

7*2221 

12.26 

7*31 

0.705 791 

7*2221 

5.41 

7*32 

0.144 040 

7*2222 

218.31 

7*32 

0.081 788 

7*2222 

286.31 

F33 

0.751 132 

7*3111 

- 6.78 

7*33 

0.758 736 

7*3111 

- 1.53 

Fm 

- 69.3948 

7*3211 

- 0.69 

7*m 

- 82.3006 

7*3211 

- 0.64 

F211 

- 3.4084 

7*3221 

0.60 

7*211 

- 3.2506 

7*3221 

0.82 

F221 

2.9664 

7*3222 

- 2.03 

7*221 

1.6407 

7*3222 

- 1.17 

F222 

- 42.4825 

7*3311 

- 2.35 

7*222 

- 47.3035 

7*3311 

2.92 

F311 

- 2.9653 

7*3321 

- 0.75 

7*311 

- 1.4453 

7*3321 

- 0.82 

F321 

- 0.2343 

7*3322 

0.47 

7*321 

- 0.2912 

7*3322 

0.65 

F322 

- 0.3804 

7*3331 

0.58 

7*322 

- 0.3408 

7*3331 

0.37 

F33 1 

- 0.6733 

7*3332 

0.50 

7*331 

- 0.6832 

7*3332 

1.37 

F332 

0.0139 

7*3333 

- 0.61 

7*332 

- 0.1290 

7*3333 

- 1.29 

F333 

- 0.7522 



7*333 

- 1.0334 
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Table 3: The CcCR Vibration-Rotation Interaction Constants, Quartic and Sex- 
tic Centrifugal Distortion Constants, and S Reduced Hamiltonian Terms for 1 
3 A' HNC and HOC+. 



Vib-Rot Constants (MHz) 
mode a A a B oP 


Distortion Constants 
(MHz) 

(Hz) 


Watson S Reduction 
(MHz) 

(Hz) 

HNC 

i 

31523.1 

-115.7 

20.0 

7 -/ 

' aaaa 

-1080.703 

10“ 3 

509.210 

Dj 

0.139 

Hj 

-0.174 


2 

2944.3 

519.5 

433.7 

T > 

'bbbb 

-0.652 

‘1* (>(,(, 

0.056 

Djk 

5.695 

Hjk 

57.640 


3 

-38642.0 

-191.0 

159.6 

7 _/ 

' cccc 

-0.494 

rc.C 

-0.124 

D k 

264.342 

10- 3 H K j 

5.567 






^ aabb 

-21.472 

10“ 3 4Va6 

15.462 

di 

-0.010 

10- 3 H k 

503.586 






T f 

‘ aacc 

-2.372 

^abb 

103.723 

d2 

-0.002 

h\ 

0.025 






^~bbcc 

-0.555 

10“ 3 <S>aac 

-9.854 



h2 

0.070 








®bbc 

-0.417 



h3 

0.020 








®acc 

5.539 












^* bcc 

-0.214 












^ abc 

120.024 





DNC 

1 

12493.5 

13.2 

120.0 

T 1 

aaaa 

-407.104 

10" 3 $0,0.0. 

122.646 

Dj 

0.113 

Hj 

-0.054 


2 

1587.0 

446.8 

368.8 

T < 

1 bbbb 

-0.610 

'I ’bbb 

0.534 

Djk 

2.561 

Hjk 

41.101 


3 

-17726.9 

-195.0 

134.1 

T f 

cccc 

-0.352 

^ccc 

-0.008 

D k 

99.102 

10- 3 H K j 

60.494 






^ aabb 

-10.855 

10“ 3 Qaab 

4.293 

di 

-0.016 

10- 3 H k 

122.545 






T f 

aacc 

-0.223 

$abb 

71.536 

^2 

-0.004 

hi 

0.089 






^ ~bbcc 

-0.435 

10“ 3 4>aoc 

-4.208 



h2 

0.159 








®bbc 

-0.578 



h3 

0.047 








$acc 

8.667 












^* bcc 

-0.078 












*^* abc 

91.369 





HOC+ 

1 

31610.2 

-11.5 

82.8 

-7 -/ 

' aaaa 

-1601.803 

10“ e <s>aaa 

1.007 

Dj 

0.127 

Hj 

-0.082 


2 

6108.6 

457.5 

403.1 

T ' 

1 bbbb 

-0.584 

'I ’bbb 

0.109 

Djk 

6.140 

Hjk 

20.640 


3 

-55101.4 

-135.0 

165.2 

T f 

' cccc 

-0.460 

4? CCC 

-0.086 

D k 

394.189 

10- 3 Hkj 

8.999 






T ' 

aabb 

-22.260 

10“ 3 3>aab 

21.388 

d\ 

-0.008 

10- 3 H k 

997.663 






T f 

' aacc 

-3.278 

&abb 

73.068 

d2 

-0.002 

hi 

0.034 






T > 

T bbcc 

-0.507 

10 -3 4>ooc 

-12.422 



h2 

0.047 








4? (,6c 

-0.241 



h3 

0.015 








4>occ 

-4.631 












^bcc 

-0.112 












4 > a6c 

96.065 





DOC+ 

1 

12393.9 

82.5 

153.0 

T 1 

‘ aaaa 

-585.740 

10“ 3 

228.978 

Dj 

0.103 

Hj 

0.031 


2 

3840.9 

386.4 

343.7 

r ' 

bbbb 

-0.543 

4*666 

0.484 

Djk 

2.728 

Hjk 

13.591 


3 

-25342.2 

-134.4 

147.4 

T f 

' cccc 

-0.331 

4* ccc 

0.014 

D k 

143.603 

Hkj 

650.885 






aabb 

-11.150 

10“ 3 4> aa 6 

6.076 

di 

-0.013 

10- 3 H k 

228.313 






T f 

‘ aacc 

-0.533 

4*a66 

47.679 

d2 

-0.003 

hi 

0.081 






^~bbcc 

-0.397 

10“ 3 4> aac 

-5.458 



h2 

0.109 








4* 66c 

-0.374 



h3 

0.037 








4>occ 

4.279 












4*6cc 

0.012 












4*a6c 

77.466 
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Table 4: The CcCR QFF harmonic and anharmonic (from VPT and VCI) fun- 
damental vibrational frequencies (in cm' 1 ) for 1 3 H' HCN, HCO + , and their 
respective deuterated forms. 


Molecule 

Description 

Mode 

Harmonic 

Previous" 

VPT 

VCI 

HNC 

a’ N— H stretch 

V\ 

3189.3 

3103 

2923. 9 b 

2922.9 


a’ N— C stretch 


1568.9 

1578 

1534.0 

1534.0 


a’ H— N— C bend 

Ds 

1165.7 

1200 

1137.1 

1134.9 




6378.6 


5591.4 

5592.9 



2u 2 

3137.8 


3037.3 

3038.4 



2u 3 

2331.4 


2253.7 

2245.0 



12 1 + V2 

4758.2 


4456.2 

4447.1 



V\ + 123 

4355.0 


4050.2 

4044.8 



"2 + V 3 

2734.6 


2666.5 

2664.3 

DNC 

a’ N— D stretch 


2329.4 


2183.5“ 

2184.8 


a’ N— C stretch 

V2 

1555.4 


1523.5 

1523.5 


a’ D— N— C bend 


891.5 


875.0 

874.1 



2v\ 

4658.8 


4232.0 

4232.3 



2 u 2 

3110.9 


3017.7 

3017.3 



2v 3 

1783.0 


1739.1 

1735.6 



12 \ + V2 

3884.9 


3710.0 

3705.8 



V 1 + V 3 

3220.9 


3055.2 

3051.7 



V2 + V 3 

2446.9 


2395.9 

2394.4 

HOC+ 

a’ 0— H stretch 

V\ 

3370.7 


3146.6" 

3146.0 


a’ 0— C stretch 

V2 

1661.9 


1630.1 

1630.2 


a’ H— 0— C bend 

v 3 

1208.3 


1168.4 

1166.8 



2v\ 

6741.4 


6080.4 

6080.4 



2v 2 

3323.8 


3230.1 

3230.6 



2v 3 

2416.6 


2306.1 

2299.5 



V 1 + ^2 

5032.6 


4774.1 

4770.6 



V 1 + 123 

4579.0 


4297.4 

4292.5 



V2 + 123 

2870.2 


2798.4 

2797.4 

DOC+ 

a’ 0— D stretch 

12\ 

2456.1 


2333.7“ 

2334.3 


a’ 0— C stretch 

122 

1643.9 


1612.8 

1612.9 


a’ D— 0— C bend 

123 

922.7 


900.9 

900.1 



2v\ 

4912.1 


4559.5 

4557.1 



2 u 2 

3287.7 


3194.8 

3194.8 



2 123 

1845.4 


1784.2 

1780.9 



12 \ + u 2 

4099.9 


3948.5 

3943.4 



12\ + 123 

3378.8 


3228.4 

3222.1 



122 + 123 

2566.6 


2517.5 

2516.3 


“Harmonic 10/(8, 2) CASPT2 frequencies from Rajendra and Chandra (Ref. 17). 
"These modes are affected by a 2v 2 = Fermi Resonance. 

“These modes are affected by a 123 + v 2 = 12 \ Fermi Resonance. 
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